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Sonar  Processing  Architectures 
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History 


Modeling  multibeam  intensity  as  a  histogram 

•  Perlovsky  (c.  1991),  Luginbuhl  (c.  1999) 

-  Interpreted  cell-level  sensor  data  amplitudes  as  histogram  counts 


•  Streit  (c.  2000),  Streit  (c.  2001) 

-  Treated  broadband  intensity  as  a  histogram 

-  Modeled  the  superposition  of  energy  from  multiple  targets  using  a 
mixture  density 

-  Extended  histogram  interpretation  to  frequency-azimuth  domain 

Direct  energy  superposition  model 

•  Ristic,  Farina,  Hernandez  (c.  2004) 

-  Used  a  model  of  the  sensor  “point-spread  function”  to  describe  the 
distribution  of  energy  across  cells  for  tracking  on  image  data 

-  Applied  a  simple  energy  superposition  model  for  developing  a  CRLB 

-  No  longer  treating  energy  distribution  as  a  pdf 


Basic  Model 


•The  basic  superposition  model 

k 

Zt={zt,l’Ztj2,--,Zt,n}T  =  Ctln  +  Yj-(Xt^  +  \ 

• _ ^ 

•The  augmented  state  J~ 

Xt  =  {C,,(x^r,(x,2)J,...,(xf  f  f 


*/  =  {  P/' » P  i  ’ It  >  Y /  }r 

•The  target  viewed  through  the  sensor  point 
spread  function 

h( x{ )  =  {h1  (x/ ), h2 (x/ ), /f3 (x/ hn (x/ )  }r 
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Estimation  Algorithm 


•  Non-Gaussian  noise 

-  Not  a  problem  for  filtering,  optimality  sacrificed 

-  Exponentially  distributed  frequency  cells  yield  Gamma  distributed 
broadband  intensities,  closely  approximated  by  Gaussian 

•  Applied  straightforward  Kalman  filter 

-  Could  use  smoother,  MLE  or  other 

-  Relatively  high  dimensionality  compared  to  traditional  trackers 

•  n-vector  measurement 

*  km+1  vector  state 

•  Covariance  decoupling 

-  If  prior  covariance  is  decoupled,  so  is  much  of  the  processing 

-  Kalman  gain  can  be  performed  with  a  km+1  vs.  n  dimensional 
inversion 

-  Output  covariance  is  fully  coupled,  but  little  performance  penalty 
seen  from  extracting  target  blocks  to  form  a  decoupled  prior  for 
the  next  update 
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Single  Target  Tracking  Results 

Fixed  Amplitude,  Varying  width  Targets 
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Improving  SNR  &  Separability 
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Non  -  parametric  model  of  target 
spectral  characteristics : 

St  ={st,i,St,2,-st,k}T 
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Outer  product  forms  model  of 
frequency  -  azimuth  image : 

FRAZ  =  hi  {xj  )T  S{ 


Parametric  model  of  spatial  location : 

•  • 

fy(x/ )  =  //  exp< 

HP/-P/)2 

2  Y t 

K,  J 

> 

•  Hold  S'  fixed,  estimate  x' 

•  Given  estimate  of  estimate  S' 
as  a  weighted  average  over  beams, 
weighting  based  on  h(xJ) 
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Summary 


•  Initialization  requires  detection,  but  tracking  does  not 

•  Superposition  model  results  in  an  implicitly 
multitarget  algorithm,  no  combinatorial  problems 

•  Simple  model  admits  simple  processing 

•  Filter  dimensionality  is  not  a  problem,  simplifying 
approximations  can  make  processing  even  simpler 

•  Provides  reliable  track  bearing  quality  outputs 

•  “Self  tuning” 

•  Tracks  over-resolved  targets  without  modification 


